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Abstract. We consider measures supported on the bi-circle and review the 
recurrence relations satisfied by the orthogonal polynomials associated with 
these measures constructed using the lexicographical or reverse lexicographical 
ordering. New relations are derived among these recurrence coefficients. We 
extend the results of [H] on a parameterization for Bernstein-Szego measures 
supported on the bi-circle. 



1. Introduction 

In this paper we continue the investigation begun in |8] on the orthogonal poly- 
nomials associated with measures supported on the bi-circle. In more than one 
variable an important consideration is which ordering to use. The usual ordering 
is the one suggested by Jackson [10] which is the the total degree ordering. This 
is natural since the addition of any new polynomials does not alter the previous 
orthogonal polynomials already constructed. However in their solution of the two- 
variable Fejer-Reisz problem Geronimo and Woerdeman [7] were led to consider 
polynomials obtained using the lexicographical or reverse lexicographical ordering 
(for an alternative viewpoint see Knese [H]). Orthogonal polynomials obtained 
using these orderings were first studied by Delsarte et al [5] who used them to solve 
the half-plane least squares problem [3] . Important in their work and later empha- 
sized in [8] is the fact that in these orderings the moment matrices have a doubly 
Toeplitz structure. This allows a connection between the polynomials obtained us- 
ing the above orderings and matrix orthogonal polynomials on the unit circle [2] 
(see also [8]). The various convergence properties of these polynomials in a strip as 
well as their connection to generalized Schur representations and Adamjan, Arov, 
and Krein theory were developed in [5], and [6]. 

Given a positive Borel probability measure a supported on the unit circle with 
an infinite number of points of increase let {4>n}n>o be the sequence of polynomials 
of exact degree n in z = e lS with positive leading coefficient having the property, 



^(e**)0 fc (e ifl )d<7(0) =hk 

T 

where Sj^k is the Kronecker delta. These polynomials are known to satisfy the 
recurrence formula [9], [12], [13] . 

4> n {z) = a n (z(j) n -i(z) - a n (/) n -i(z)), n > 1 

where <j> n {z) — z n <ft n (l/z) is called the reverse polynomial. The a n are called the 
recurrence coefficients and a n — ? -1 where k n is the leading coefficient of <p n . 
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From the orthogonality properties of <f> n and <fi n -i it is not difficult to obtain the 
relation 

f = al(l-\a n \ 2 ). 

Thus \a n \ < 1 and given a n , a n can be computed. The recurrence coefficients play 
an important role in theory of orthogonal polynomials on the unit circle as can be 
seen by Verbhmsky's Theorem [12] , 

Theorem 1.1. Let a be a Borel probability supported on the unit circle with an in- 
finite support, then associated with a is a unique sequence of recurrence coefficients 
{ a n}^Li with \a n \ < I. This correspondence is one-to-one. 

A useful characterization theorem is the following, 

Theorem 1.2. Let a be a Borel measure supported on the unit circle. Then a is 
absolutely continuous with respect to Lebesgue measure with density ^ where 
p n (z) is a polynomial of exact degree n in z with *p n (z) nonzero for \z\ < 1 if and 
only if cti = 0, for i > n. 

Measures of the form given in the above Theorem have come to be called 
Bernstein- Szego measures [P2"] . 

In [8 a parameterization of the two variable trigonometric moment problem was 
introduced in an attempt to be able to find extensions of the above two Theorems 
to the two variable case. Here we continue the study of the algebraic properties of 
this problem. In section 2 we call together the results needed. In particular begin- 
ning with Borel measures supported on the bi-circle whose moment matrices are 
positive we construct orthogonal polynomials using the lexicographical or reverse 
lexicographical ordering. Then recurrence relations satisfied by these polynomials 
are displayed and some properties of the recurrence coefficients are noted. Here 
the parameterization discussed above is introduced and a two variable analog of 
Verblunsky's Theorem is presented. In section 3 we develop new equations between 
the recurrence coefficients that shed light on how these coefficients are related to 
each other. Some of these recurrence relations are used in section 4 to develop an 
algorithm different from that given in [8] which allows us to make more precise the 
construction of the parameters left undetermined in Theorem 7.9 of [8]. 

2. Preliminaries 

In this section we collect some results that will be used later. As noted above 
we will use the lexicographical ordering which is defined by 

(k,£) <io X (h,£i) ^> k < ki or (k — ki and t < £i), 

and the reverse lexicographical ordering, defined by 

(M) < rev lex 

Both of these ordcrings are linear orders, and in addition they satisfy 

(fc, £) < (m, n) (k + p, £ + q) < (m + p.n + q). 

In such a case, one may associate a half-space with the ordering which is defined by 
{(k, I) : (0, 0) < (fc, I)}. In the case of the lexicographical ordering we shall denote 
the associated half-space by H and refer to it as the standard half-space. In the case 
of the reverse lexicographical ordering we shall denote the associated half-space by 
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H. Let a be a positive Borel measure support on the bi-circle z — e ie , w = e' 1 ^ 
with Fourier coefficients, 

Cfcj = / e- lke e- l ^da{0^). 
Jt 

We now form the (n + l)(m + 1) x (n + f)(m + f) moment matrix C n ^ m using 
the lexicographical ordering. As noted in the introduction it has a special block 
Toeplitz form 



(2.1) 



Co C_i 
Ci Cq 



c_ n 

C_„+i 



_C n C n -i ■ ■ ■ Co . 
where each C is an (m + 1) x (m + 1) Toeplitz matrix as follows: 



(2.2) 



Ci = 



Ci,0 



i = —n, . . . , n. 



Thus C„ jm has a doubly Toeplitz structure. If the reverse lexicographical ordering 
is used in place of the lexicographical ordering, we obtain another moment matrix 
C n , m where the roles of n and m are interchanged. Throughout the rest of the 
paper we will assume that C„ jTrl is positive definite for all < n, m. 

We now compute orthogonal polynomials associated with a. We begin by order- 
ing the monomials z % w\ < i < n, < j < m lexicographically then performing 
the Gram-Schmidt procedure using this ordering. Define the orthonormal polyno- 
mials <fi l n m (z,w), < n, < m, < / < m, by the equations 

(2.3) 



/ 

J 

JT 



x z l w 3 da = 0, < i < n and < j < m or i — n and < j < I, 



n,m*rn,m 



and 
(2.4) 



n, m 



(z,w) = K'LiZ n ™ l + 



E 

(»>j)<lex(ra,() 



With the convention k nm x > 0, the above equations uniquely specify 4> L n m . Polyno- 
mials orthonormal with respect to a but using the reverse lexicographical ordering 
will be denoted by (f> l n m . They are uniquely determined by the above relations with 
the roles of n and m interchanged. 
Set 







' 6 m 




- z n w m - 


(2.5) 


<f>n,m = 


±m-l 
_ 4*n,m _ 




z n w m-l 
1 
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where the (m + 1) X (n+ l)(m + 1) matrix K n ,m is given by 



(2.6) 



K„ 



h,n,m 
n,m,m 



h.n,m—l 
" J n,m,m 
i n,m— 1 



As indicated above denote 



(2.7) 





h n,0 


7tt-l,m 
K n,m,0 








- w m z n - 




— K 
"71,771 


w m z n-l 


6° 




1 









1,0,0 
71,771,711 
r 0,0 



,0,0 

€ ra,m,0 



where the (n + 1) x (n + l)(m + 1) matrix K n ^ m is given similarly to (|2.6p with the 
roles of n and m interchanged. For the bivariate polynomials 4> l n m ( z i w ) above we 
define the reverse polynomials 4> l nm (z, w) by the relation 



(2.8) 



0' m (z,w)=z n w m 4>L m (l/z,l/w). 



«-r, 



With this definition cf> n m (z,w) is again a polynomial in z and iu, and furthermore 

T 



(2.9) 



m— 1 

71,711 





71,711 



An analogous procedure is used to define 

To find recurrence formulas for the vector polynomials <& n ,rm we introduce the 
notation for every vector valued polynomials X and Y, 



(2.10) 



(X,Y)= / X(z,w)Yi(z,w)da, \z\ = l = 
Jt 



The following recurrence formulas which follow from the orthogonality relations 
satisfied by <fr n ,m and $„. m were proved in [8]. 

Theorem 2.1. Given {$ ra ,m} and {& n ,m}, < n 7 < m, the following recurrence 
formulas hold: 



(2.11) 
(2.12) 
(2.13) 

(2.14) 
(2.15) 
(2.16) 



i rp 

^71,771^71,711 — Z*&n— l,m E n m < E > n _i_ m , 

"^n,™ + ^n,m^n,m(^-Ti,m) ^71,771 = ^n,m z< ^n-l,"i! 

1^71,771^*71,771 ^71,771—1 ^71,711^71—1,711? 

^71,771 -^71,771^71,771 ~t~ 771 ^ 71 , 771 — 1 1 

^71,771 — ± n,m^n,m < V A 71,771/ ^ 7i,m— ll 
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where 






(2.17) 






(2.18) 






(2.19) 




= (*»,*,-!, *»-l,m) GM m ' n , 


(2.20) 


V 

- 1 - n.m 


= (*»,„,_!, $„,„,) £M m ' m +\ 


(2.21) 


AT 1 

n,m 




(2.22) 


r 1 

n,m 


= (W*n, m - U * n , m ) eM m ' m+ \ 


(2.23) 




— /(b <T> \ c )ifm+l,n+l 


(2.24) 


i 1 

n.m 


— (% T <J> \ £= JVf m+1 ' n+1 



Here M 1 ^ denotes the set of i x j matrices with complex entries. Equation (12. lip 
was first found by Delsarte et. al. [2]. 

Remark 2.2. Formulas similar to (|2.11[) - (|2.16p hold for $> n>m and will be denoted 
by ( 12.lip -f 12.16p . Throughout the rest of the paper we use the same notation to 
denote the tilde analogues of existing formulas stated for § n ,m- 

Examination of equation (|2.13p shows that the entries of T n ^ m are zero 

for i > j with the (i, i + 1) entries positive. Likewise equation (12. 14|) implies that 
entries of m are zero for i > j with the (i, i) entries positive. 

From the definitions of /C, JC 1 , I, I 1 and their tilde analogs it is not difficult to 
see the following relations, 

(2.25) K"n,m — K„ m , In,m — In,mi 

(2.26) In,m = {^n,m) ) ^n,m Q^"n,m) • 

Also the recurrence relations yield, 

(2.27) A n ^ m A^ n ^ m = I m — E n . m E^ ml 

(2.28) ^ n,mX\i^ m ^"n,mK"n,mi 

(2-29) r„ m (r„ jm )^ = i m — ic n ^ m {fc n , m )\ 

(2.30) In,m,In : m ^n,mJ^n,m = Im+lj 

(2-31) In,m(In,m) (J-n,m) ^n,m = 

The definition of $> n , m implies that A n ^ m is an upper triangular matrix with 
positive diagonal entries. Thus it may be computed from E n m using a Cholesky 
decomposition of equation (I2.27p . A slightly more involved analysis [8] shows that 
r n ,m may be computed using a Cholesky decomposition of (|2.28p . 

The above recurrence formulas also give pointwise formulas for the recurrence 
coefficients. In order to obtain these formulas we define the m x m + 1 matrices 
U m and as 

(2.32) U m = [0, I m ] , 
and 

(2.33) U m = [I m , 0] , 
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where I m is the to x m identity matrix. From equations (|2.39|) and (|2.40[) we write 

$™(z) = $™^" + $I™ 
(2.34) *»(«,) 
then the following relations hold: 



$ m z n 4 
-- w n 

rn.m w 



' 71,71 — 1* 



(2.35) 
(2.36) 
(2.37) 
(2.38) 



-V T F 



-T 1 I 1 F 1 

n.m n.m ti.rai 



where F = <E>™ r/ T f$' i_1 1 -1 and F 1 = <E>™ fr/ 1 l T ('4> n_1 1 -1 

wucic J- ^m.m u n V*m,m) ' x n.m ^m.mX^n) \^m,mJ 

Equations (|2.11J and (|2.12|) are a consequence of the relation between & n ,m and 
the matrix orthogonal polynomials associated with the (m + 1) X (m + 1) matrix 
measure M m given by 

t 

dM m {6) = 





- w m- 




- w m- 










7T 


_ 1 _ 




1 



where u> = e l 'K Given E^ m equation (|2.1ip allows the computation of $i !m along 
the strip < i, < j < m. More precisely if we write 

(2.39) * n , m (z, w) = $%{z)[w m , • . • , 1] T , 

then the $™ are a sequence of matrix polynomials of degree i in z satisfying 



$T{z)dM{6){$T (2))t = I m+l 8i, 



3 ' 



where 7 m +i is the (m + 1) x (to + 1) identity matrix and (5ij is the Kronecker 
delta. Thus (|2.11[) and (|2.12[) follow the recurrence formulas satisfied by matrix 
polynomials orthogonal on the unit circle. Equation (|2.5[) implies that the coefficient 
of z l in $™, $™ , is an (to + 1) x (m + 1) upper triangular matrix with positive 
diagonal entries. Similar statements hold for 

(2.40) ®n,m{ Z M = $n MIA • • • , If- 



In contrast given ,and $ij for (z, j) = (n — 1, to) or (n, to — 1) equations (|2.13D 
and (|2 . 14[) allow the computation of $ n ,m- 

As noted in [8] there is a lot of redundancy in the coefficients of the above equa- 
tions. If we have all the Fourier coefficients in the notched rectangle {(i,j),0 < 
i < n, < j < to} \ (n, m) then the polynomials $„. m _ i, 'I'n-im can be computed. 
Notice that only two new Fourier coefficients are required to compute all the poly- 
nomial <& n ,m whereas /C„ jm and /C* m are both m x n matrices. This led in [8] to 
the introduction of parameters Uij i > 0, w-i.-j = j such that 



(2.41) 
and 

(2.42) 



-n,-m = (^n,m) 1,1 71 > °) m > °i 



(e™) T (^ 1 )- 1 /c„, m (($^)t)- i < 



( e m) ^Wi,me ra /(fc, 



i!,n 



"m.O 



n,m— 1,0 n—l,m 



), n > 0, TO > 0, 
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where is the m-dimensional vector with zeros in all its entries except the last, 
which is one, and k™' m —i and h%!\ m are the leading coefficient of 4>n.m-i and 
4>n-i m respectively. The last equality was obtained using the upper triangularity 
of and $„7^, and equations (12.51) and (|2.7|l . In terms of inner products the 

parameters can be written as 



(2.43) «_«,_,= / Cm-i (*, ™)1=U(^ «0<M*, 0), ^ = e ie , «; 
and 



(2.44) u_„, m = / < m _i(«,tu)$_i m (z, W )da(0», z = e M , w = e* 

JT 2 

Here 0° m _ 1 and m-1 have leading coefficient one. Since /C n , m and K\ m are 
contractions the parameters must satisfy the constraints 

\v>n,m\ < 1 

and 

, n,0 7m,0 I I < i 

With this the following Theorem was proved in [8] , 
Theorem 2.3. Given parameters u^j E C, < i, U-ij — construct 

• scalars and Eo,j> 

• matrices ICij, i > 0, j > 0; and 

• numbers (e[) T Hfje[, i > 0, j > 0. 

// 

(2.45) u ,o > 0, |^o| < 1, \E ,j\ < 1, H^j-ll < 1, and (H? tj )i,i < 1, 
i/ien iftere exists a unique positive measure a supported on the bi-circle such that 

(2.46) J $ iim $t m dcr = Jij-J^+x and jf ^dcr = SijI^ 

The conditions (|2.45[) are also necessary. 

The numbers (ff? -Ji^ are given by equation (5.14) in [5] . 

This Theorem is the two dimensional analog of Verblunsky's Theorem discussed 
in the introduction. A polynomial p is of degree (n, m) if 

n m 
i=0 j=0 

with p„ iTO ^ 0. 

When dcr = ^ ^ ^jp where p n ,m is a polynomial of degree (n, m) with p"n.m 

stable (i.e. Vn,m 7^ 0, |z|, < 1) then more can be said. 

Theorem 2.4. Let fi be a positive measure on the bicircle. Then ji is purely 
absolutely continuous with respect to the Lebesgue measure and dpi = i7r i® p d ^ p , 
where p n ,m is a polynomial of degree (n, m) with p n ,m stable if and only if 



i+i- 



(a) K n> j = 0, E n -x,j+x = 0, and it„j+i = 0, j > m: 

(b) Ki >m = 0, Ei, m -x = 0, and Ui, m = 0, i > n; 

(c) mj = 0, \i\ > n, j > m. 
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It was actually shown that in this case u- n j, U-i }7n , u n -i,j+i and Ui, m -i are 
equal to zero for j > to, i > n. How to compute the remaining parameters and 
how they were related to Uij, where \i\ < n and j < m, was not indicated and it 
the subject of the remaining sections. 

3. Relations for E n . m , tC„. m and K} nm 

In order to prove Theorem 12.31 it was necessary to show that most of the entries 
in K n . m and K„ m could be computed knowing the recurrence coefficients on the 
(n — 1, to) and (n, to — 1) levels. These relations will be augmented by the following 
new relations which will be used later to compute coefficients on lower levels from 
those on higher levels. 



Lemma 3.1. For n > and to > 

' r /+ ! ' "i" /v ! 1 I ! .1,1 + 



(3.1) E nm — Y n —x,m+X^n,m+x{Pn—l.m+x) + )Cn—l,m+l(J^-n—l ,ra + 1) ■ 



Also for n > and to > 0, 

(3.2) E n , m = r n+ i iTO _i£ , n+ i iTO (r^ +lim _ 1 ) T + ^ n+ i im _i(^ +ljtn _ 1 ) T . 

Proof. From the definition of E n ^ m , 

(3-3) E n>m = (z$ n -i, m ,QZ_ ljfn ) 

eliminate <5> n -i,m using equation (|2.13|) to obtain 

E n ,m = r n — ljm+l^^n-ljm+lj l,m) 

With the use of the reverse of (|2. 14p the first integral on the right hand side of the 
above equation can be rewritten as 

(z®n-l,m+l, *n-l,m) = { z ®n-l,m+l, ^n-l,m+l)(^n-l,m) T 

+ ($n-l,m+l, *n-2,m+l)(^n-l,m+l) T 
= ^n,m+l(r^-i. m +i) T - 

Equation (]2 . 1 7[) and the orthogonality of $ri-i.m+i to $„_„ im+ i has been used 
to obtain the last equality. The result now follows by taking the transpose of 
equation (|2.21j) . Equation (|3.2[) follows in a similar manner using the tilde analog 
of the above equations. □ 

Lemma 3.2. 

(3.4) Y n — i rn E n m I n _^ m — A nra —\KL nra ^Cn-l,mfn-l,m 

and 

(3-5) T^ nrn —iE n ^ 7n I nnl _-^ — A n —\ ra ]Q nvrL /Cyi,m— l-^n.fn— 1 

Proof. To obtain (|3.4[) note that equations (I2.17P and the reverse transpose of 
(1239 give, 

E n ,mln~l,m = ( z ®n-l,m, $n-l,m) ■ 

Multiplying the above equation on the left by T„_i !m then using the reverse trans- 
pose of equation (|2.13[) yields, 

— l,m) fcn— l,m 
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The second integral in the above equation evaluates to T^—i m - Substitution 
of equation (|2.11[) in the first integral to eliminate 2$ ra _i im _i then using equa- 
tion (|2.19p yields equation p.4p . The argument for equation p. 51) follows in an 
analogous manner using the tilde analog of the above equations. □ 

Finally 
Lemma 3.3. 

and 

d 7) P E ff 1 ) T - (K 1 ) T A T - r t (K 1 ) T 

V°- 1 I I ii,m-l"n,i"l 1 n,m—l) ~ V^n.mJ -"-n— l,m ^ n,m—l v^n.m — 1 ) 

Proof. To obtain (JO) use equations (l2~T7)) and ( 12~T5)) to find 

^n—l,m^n,m ~ (z&n— l,m> ^ n— l,m)" 

Multiplying the above equation on the left by the transpose of r*_ x m then using 
the reverse transpose of equation ()2.14|) yields 

- ($«-l,m, $n-2,m)(/C I 1 l _ l m ) T . 

The second integral in the above equation evaluates to f m . Substitution of the 
reverse transpose of equation (|2.11[) in the first integral then using equation ( 12.2ip 
yields equation (|3.6|) . As above equation (|3.7[) follows a similar argument using the 
tilde analogs of the above equations. □ 

With these recurrences we can prove a strengthening of Lemma 7.5 in [5] 

Lemma 3.4. If Eij = 0, then the first column offCj j is equal to zero, in particular 
Uij = 0. If Eij and ICi-i j (]€}_! j) T are zero, then so is Eij~\. IflCij, Eij-\, 
and Uij are zero, then Eij = 0. Likewise if Eij — 0, then the first row of ' Kj ■ is 
equal to zero. If Ei j, and K/ i ^_ X K\ are zero, then so is IflCi.j, Ei-ij, 

and Uij are zero, then Eij = 0. 

Proof. If Eij — 0, then (|3.6p and the triangular structure of Ti-ij show that the 
first column of K\j is zero. If Eij and /Ci_i J (/C^_ 1 -) T are equal to zero, then p.ip 
shows that E%j—\ = 0. Equations (3.50) and (3.51) in [8] are 

(^•8) ^ n — l.m-En^m -^-71,171— l^- , n,7n(.^n / — l,m) ^n,Tn— X^n — l,m3 

(3.9) En,m (r n -l,m) = In — l,m{fc n ,m) ^-n,7n-l ~^ ^n-l,m^n,m-l- 

Thus if /Ci.j and -Ej-ij are equal to zero (|3.8[) and the fact that Eij is symmetric 
shows that all its entries are equal to zero except for [Eij]^i^y Equation (13. 9[) 
and the assumption that Uij = give that this entry is equal to zero also. The 
remaining statements follow in an analogous fashion using the tilde analogs of the 
above equations. 

□ 
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For the next lemma we recast equations ()2.5|) and (|2.7|) as, 

m 

(3.10) $ n , m (z, W) = J2 Ln, m W l [z n , . . . , 1] T 

i=0 

and 

(3.11) $ n , m (z, «;) = £ L^^K", . . . , 1] 

i=0 

From this we have, 
Lemma 3.5. For n,m>0 

(3-12) I n ,m = L n ^ m (<& mm ) 

and 

where J n is the (n + 1) x (ri + 1) matrix with ones on the reverse-diagonal and zeros 
everywhere else. 

Proof. The first formula follows from equation (|3.10[) . The second equation can be 
seen from the computation, 



t n , m (z, wf = Ll m w m -'J n [z n 7 . . . , 1] 

i=0 



□ 



The above results give formulas for the parameters. From equations (|2.35p . 
(12371) . and (gIT2j) we find 
(3.14) 

where we have used the fact the Uje^ — Gj+i m tne ^ as * equation. Likewise 
equations (j2~31)l) . (|2~3"gj) . and (|5TT5|) show, 

(3-15) (/C^.)i,i = -(^ 1 ^($l l )- 1 L^J,(L^) T ((|.}- 1 )t|»}- 1 )- 1 ) M . 

4. Construction of the Parameters 

Using the results above we are now able to compute the remaining parameters 
from those given in the rectangle < i < n, < j < m. We first show that all 
the polynomials $>i t j and can be computed for i > n and < j < m. To see 
this suppose that we are given j and for < j < m and also conditions 
(a), (b), and (c) of Theorem 12.41 are satisfied. From their defining properties we 
see that 4>™ m — 0™ m . Lemma l3~4l shows that (a) and (b) imply that £^ m = for 
i > n so that from equation (|2.1ip = z l ~ n $ n ^ m and because £^ m _i = 0, i > 
n, $i, m -i = 2 4_ ™$ JliTO _i. If /Cj iTn = it follows from the triangularity of T^ m 
and equation (|2.28[) that IY,„ = U m . Likewise Ti^ m — C/j. This implies through 
equations f~ l2~T3f and f~ l2~23"|) that 

- _ Vn,m 
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for i > n which gives all of <£>i, m for i > n. Set it) = in ( 12.12[) then utilize 



4^ 

equation (|2.40[) and the fact that 3>"(0) is invertible to obtain 

(4.i) it m l i>m (ir m )-i = -^MMKM-'f = B) n . 

tr ____ 

To find Ai t , n use the orthogonality properties of $i !m and i >m in equation f * 12.12|) 
to find, 

T — R i ('RMt — Jt 4 . 

Since Ai jm is upper triangular with positive diagonal entries it may be computed 
using the lower Cholesky factorization of the left hand side of the above equation. 
Using this in ( 12.12[) allows us to compute $i jm _i for i > n. In an analogous 
fashion may be computed for i > n and < j < m — 1. Now Ti t j and T+ j, 
i>n, 0<j<m — 1 may be computed from from equations ( 12.35|) and (~ !2.36j) 
respectively. With i = m — 1 we find from (|3.1j) since E n +x m -i — that, 



which gives E n +i.m-2 because by assumption K. n +i iTn -i and /C* , 1 m _ 1 are known. 
Since A n+ i_ m ^2 may be computed from i? n+ i jm _2 using the upper Cholesky fac- 
torization of (|2.27p we obtain <& n +i.m-i from (|2.11|) . By induction we see that the 
above argument gives E n+ i^, i = 0, . . . , m — 3 from which An+i,i may be com- 
puted and then <5> n+lti . Using equations (|3.12j) . (|3.13|1 . (|2.35|) - (|2.38j) allows us to 
compute r„ +M ,ri +M , I n+ i ti , 1^+1,%, fcn+i,i, and K} n+li for < i < m - 1. With 
equation (|3.ip and the coefficients just computed we repeat the above argument for 
level (n + 2, i) and by induction i > n, < j < to — 1. 

We summarize this with 

Lemma 4.1. Given (a), (b), and (c) of Theorem \'2A\ as well as <& n ,j an d ^n,j f 0T 

< j < m then and $i j for i > n, < j < m can be computed recursively. If 
$i iTO and $i !?rl are given then and for < i < n, j > m can be computed 
recursively. 

We now use the formulas (|3.14l) and (|3.15l) which give for i > n and 

1 < j < m and Uij for i > n and 1 < j < m — 1. For j = we have from [8] that 

(4-2) = 

which gives the parameters in the strip i > n, < j < m — 1. To compute 
the parameters in the strip j > m, < i < n — 1 equations (|2.25[) and (I2.26|) 
show that we need only interchange i with j and the matrices associated with 
the lexicographical ordering with those associated with the reverse lexicographical 
ordering in equations (|3.14[) . (|3. 15|) . and (|4.2[) . This leads to, 

Theorem 4.2. Suppose a is a positive Borel measure supported on the bi-circle 
with parameters Uij . If 

(a) Ui t -j, 1 < i,j < n and Ujj, < i,j < n give /C„. m = 0; 

(b) ttjj = /or i = n — 1, j > m, i > n, j = m — 1, \i\ > n, j — m, \i\ > n,j > to; 

(c) for i > n, < j < m — 2, Uij are equal to the left hand sides of jff.J^| ). 
H3.15\) . or \4-<fy 'computed using the above algorithm; 
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(d) for j > m, < i < n — 2, m ,j are egwaZ to i/ie left hand sides of ( ]3.14\ ). 
( ]3.15\) . of ( ]4-2ty computed using the above algorithm, 

then a is absolutely continuous with respect to Lebesgue measure with density 1 p 
where Pn^ m is of degree (n,m) with p n ,m{z,w) stable. p n ,m is unique up to multi- 
plication by a complex number of modulus one. 

Proof. If ICn^m = then from Theorem (7.3) of [8] or Theorem (10.1) of [11] there 
exits a polynomial p n ^ m of degree (n, m) with Vn.m stable such that the measure 
dp = A ^2\p^ 1 2 has parameters given in (a). By Theorem (I2.4[) and the algorithm 
above we that this measure has the same paramters as in (b) and (c). Thus by 
Theorem §3J da = A J\v* m? - 

□ 
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